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We study the thermodynamical properties of the quantum 
kicked rotator, coarsened by an external fluctuation with a 
weak intensity D, by means of the Tsallis entropy with a 
changing entropic index q. The genuine entropic index, cor- 
responding to given values of D and h is that making the 
Tsallis entropy increase linearly in time, and it is proved to 
become q < 1 for suitably large values of h: This indicates 
a subdiffusional regime which, in turn, signals the occurrence 
of quantum localization. Thus the process of Anderson local- 
ization is shown to be compatible with a thermodynamical 
representation provided that a non-extensive form of entropy 
is used. 
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The study of the connection between dynamics and 
thermodynamics is making significant progresses along 
the lines of the seminal work of Krylov Interest- 
ing attempts are currently being made at relating the 
Kolmogorov-Sinai (KS) entropy Q to the thermodynam- 
ical entropy. Of remarkable interest are the work of Gas- 
pard, relating the KS entropy for a dilute gas to the 
standard thermodynamical entropy per unit volume of 
an ideal gas j^], and the more recent paper by Dzugutov, 
Aurell and Vulpiani who express the KS entropy of 
a simple liquid in terms of the excess entropy, namely, 
the difference between the thermodynamical entropy and 
that of the ideal gas at the same thermodynamical state. 

On the other hand, the realization of Boltzmann's 
dream H implies the Gibbs entropy 



s G (t) 



dXp(X,t)log(p(X,t)) 



(1) 



corresponding to a given nonequilibrium initial condition, 
to quickly increase towards the equilibrium condition es- 
tablished by Boltzmann's entropy. The time evolution 
of the Liouville density p(t) is, on the contrary, unitary, 
thereby making the Gibbs entropy time independent. As 
a consequence, in accordance with the view of Krylov @] 
a coarsening process is required. Since, the second princi- 
ple n requires the system to be thermally isolated from 
the environment, and not necessarily totally isolated, we 
deem to be acceptable to perturb its quantum equation 
with a stochastic force of weak intensity D [Q: With 
heuristic arguments these authors prove that this 
weak perturbation makes the Gibbs entropy of Eq. (jl]) 
reach the maximum value established by Boltzmann's 



entropy in the short time scale 1/A, with A denoting the 
Lyapunov coefficient. This heuristic prediction has been 
more recently checked numerically by Pattanayak and 
Brumer |8[ who used rather than the conventional Gibbs 
entropy another form of extensive entropy, the Renyi en- 
tropy |J. Thus, the entropy increase has a rate fixed by 
A2, where A2 is a generalized Lyapunov coefficient coin- 
ciding with A in the ideal case of strong chaos. 

The important problem of how to deal with the case of 
non extensive thermodynamics has been addressed about 
ten years ago by Tsallis fl0|| , with a generalized entropy 
that, expressed in terms of the Liouville density pit), 
reads: 



1 - / dXp(X,t) 



(2) 



As an important benefit this entropy makes it possible to 
extend the KS method of numerical analysis to the frac- 
tal dynamics [^TJ.The authors of Ref. jllj generalized the 
KS entropy by replacing the conventional Shannon ex- 
pression, formally equivalent to the structure proposed 
by Gibbs, with the non-extensive form of Tsallis. In the 
special case q = 1, when Eq.((|) yields Eq.(@), 

the gener- 
alized KS entropy, denoted by K q , becomes identical to 
the conventional KS entropy ||,|| . 

We establish a connection between K q and the gener- 
alized entropy of Eq. (^|) with the same coarsening pro- 
cedure as that originally proposed by Zurek and Paz M . 
However, to adapt it to the special case of fractal dynam- 
ics we found to be convenient to use the intuitive picture 
by Zaslavsky pjj . We adopt an initial condition compat- 
ible with classical physics so as to generate a trajectory 
that will explore in due time a given region of the phase 
space according to the prescriptions of classical mechan- 
ics. Let us call the portion of this region not yet 
visited at the given time t by the classical trajectory. It 
is evident that in general this quantity will be a decreas- 
ing function of time. For generality, on the basis of the 
arguments used in Ref. [jll| to generalize the KS entropy, 
we express the general form of this time evolution as: 



R(t) = 



R(0) 



[l + X q (l-q)t]T^ 



(3) 



where X q denotes a generalized Lyapunov coefficient. Let 
us imagine that at the very same moment when the ini- 
tial conditions are established a weak fluctuation process 
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of intensity D is switched on. In the absence of fluctua- 
tions a set of trajectories, with the same initial condition, 
will behave as a single trajectory, even if this single tra- 
jectory will be made to look erratic-like by deterministic 
chaos. The action of an even extremely weak fluctuation, 
in the presence of chaos, will force the trajectories of this 
set to spread over the proper phase-space region even if 
all of them are assumed to be given exactly the same 
initial condition. This weak stochastic force serves the 
purpose of defining in a precise way the time at which 
the exploration of the phase space can be regarded as 
being completed. We call this time tcG and we define it 
by means of the following relation: 



R(t CG ) = Dt CG - 



(4) 



We note that the solution of this fundamental equation 
depends on the nature of the dynamical process under 
study. In the case of strong chaos JllJ we have q = 1 
and, consequently, 



tcG 



-log(l/V£>), 
M 



(5) 



which coincides with the estimate made by Zurek and Paz 
to determine the onset time of thermodynamics. How- 
ever, in the case of highly correlated dynamical processes, 
with q < 1 as in the case here under study, we obtain from 

Eq. m 



where 



1 

tcG « -fig, 



9 = (l-q)/(2-q). 



(6) 



(7) 



This result gives the misleading impression that the 
thermodynamical regime is dramatically postponed in 
the case of dynamical processes with no time scale. It is 
important to stress that this postponment only concerns 
the onset of ordinary thermodynamics. In fact, a weak 
stochastic force perturbing the dynamics with long-range 
correlations responsible for anomalous diffusion has the 
key effect, as shown in Ref. of producing a crossover 
from anomalous to ordinary diffusion at the time scale 
of Eq.(|6|). However, according to the perspective estab- 
lished by Tsallis jlOj , also the regime of anomalous diffu- 
sion is given a thermodynamical significance if the proper 
cntropic index q =/= 1 is used. Therefore we expect that 
the coarsening stochastic force producing the crossover 
from anomalous to ordinary diffusion at the time scale 
of Eq . (^) results in the fast increase of the entropy S q of 
Eq.(|2). This means that a transition from the regime of 
Eq. (|) to that of Eq. (g) (with q ^ 1) occurs and con- 
sequently that the entropy increases with the rate 1/X q . 

In this letter we support this conjecture in a case 
where the origin of long-range correlation is quantum. 
As pointed out by Zurek and Paz [Q, the time scale of 
Eq.(S) has to be compared to 



(8) 



which corresponds to the onset of long-range quantum 
correlations. The transition from quantum to classical 
physics takes place when \TB is made larger than h. 
Rather than focusing our attention on values of h so small 
as to ensure the classicality condition || , we explore the 
transition region, and the long-range correlations of the 
region h > \/D.The case study is given by the quantum 
kicked rotor. This means the Hamiltonian 



H(t) 



21 d9 2 



e cos 9 5{t-nT), 



(9) 



where I is the rotor moment of inertia and 9 is the ro- 
tor angular coordinate. This Hamiltonian results in the 
quantum map 

UT d 2 

{6\<p{n + 1)) = cxp(z — — )exp(-ikcos9)(6\(l>(n)), 

(10) 

implying the classical control parameter K to be ex- 
pressed as: K = Yfk. 

Here we adopt two different numerical methods to eval- 
uate the entropy time evolution. The former is the same 
as that designed by Pattanayak and Brumer, based on a 
Fourier transform technique ||. The latter, henceforth 
referred to as a statistical ensemble (SE) method, is based 
on the study of the time evolution of N independent sys- 
tems, all of them driven by Eq.(10), whith a random in- 



gredient artificially added by replacing exp y—i^jf-^j with 

exp {^ ^ P+ 2ih^ )' wri ere F(n) is a computer generated 
random process. This computed generated random pro- 
cess is introduced to play the same coarsening role as the 
earlier mentioned stochastic force of the picture of Zurek 
and Paz Q| . For this reason the variance of this random 
process F(n) is called D, is assigned very small values 
and is referred to throughout as intensity of the stochas- 
tic force. The statistical density matrix is then evaluated 
by means of the ordinary statistical average: 



p(n) 



1 N 

-Y\Mn))(Mn)\, 



(11) 



usually setting N — 10. The same method has been used 
in an earlier paper Jig] and proved to produce results 
equivalent to those of the literature on this subject. It 
is important to stress that the ingenious Fourier method 
of Ref. H establishes a connection with the system dy- 
namics unperturbed by the stochastic force, much in the 
spirit of the linear response theories. The method here 
adopted, on the contrary, allows us to go beyond the 
second-order approximation, and makes it possible to de- 
tect, as we shall see, the saturation effects determined by 
the steady action of the stochastic force. 
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FIG. 1. The Renyi entropy as a function of time. The 
values of Ti adopted are from the the top to the bottom curve: 
h = 0.01, 0, 05, 0.1. We also set K = 7.1 

In Fig|l| we show the increase of Renyi entropy as a 
function of time for different values of h. This calcula- 
tion was done using the same method as that of Rcf. ||] 
with a choice of parameters corresponding to the quan- 
tum condition that we plan to explore by using the SE 
method. We checked that the SE method yields tha same 
results as those illustrated in Fig.]!], if the same parame- 
ters are used and we set \[~D — 0.002. We note that Pat- 
tanayak and Brumer Q considered also values of Ti much 
smaller than the value h = 0.01 here considered. This is 
the reason why we depart from the regime of exponential 
increase, searched and found by those authors, and we 
rather get the power law increase: S(t) — cons t a with 
a 2.5. This property is compatible with the numerical 
results obtained by Pattanayak and Brumer at their large 
values of h. These authors, however, seemingly did not 
consider this behavior to be thermodynamically relevant. 
On the basis of our earlier remark, on the contrary, we 
are here in a position to disclose the thermodynamical 
nature of this regime, in spite of the fact that one might 
judge it dominated by the quantum mechanical coher- 
ence, and by the Anderson quantum localization |^] as 
well. 

We are convinced that at h = 0.01 we are in the pres- 
ence of a condition equivalent to that of fractal dynamics 
with no time scale, and thus corresponding to an entropic 
index q =/= 1. However, with increasing the intensity of the 
stochastic force we can provoke a transition from quan- 
tum to classical physics which, as earlier pointed out, 
must take place at \f~D = ti. Fig|| shows the following 
very remarkable property. At small values of D the curve 
Sc(t) is characterized by a positive second-order deriva- 
tive with respect to time (convex curve) . With increasing 
the noise intensity a transition from concavity to convex- 
ity is produced at about the value of D (\f~D w 0.01) 
corresponding to the passage from quantum to classical 
physics. We assume that the transition from convex- 
ity to concavity has a thermodynamical relevance, and 
this leads us to find still more interesting results by plot- 



ting the Tsallis entropy as a function of time for different 
values of the entropic index q. This is shown in FigJ|, 
which illustrates the remarkable property that decreas- 
ing q has effects similar to those obtained by increasing 
D (see Fig||). The critical value of q corresponding to 
this transition is q « 0.33. 



FIG. 2. The Gibbs entropy as a function of time for 
h — 0.01, K = 7.1 and different values of the noise inten- 
sity D. The values of D changes from the bottom to the top 
curve as follows: \JT5 = 0.002, 0.006, 0.01, 0.02, 0.05. 




FIG. 3. The Tsallis entropy as a function of time. The 
values of the parameters used are h = 0.01, K — 7.1 and 
\f~D = 0.002. The entropic index q changes from the top 
to the bottom curve as follows: q = 0.155,0.205,0.255, 
0.280,0.305,0.330,0.335,0.405, 0.530,0.555 and 0.805. 

To make our finding still more impressive we plot in 
Fig.|| the saturation time t$ as a function of the noise 
intensity D. To derive the results of Fig.[| we defined 
the saturation time as the time at which the curve Sc(t) 
reaches the 50% of its maximum value. The remarkable 
fact illustrated by this figure is that the full line corre- 
sponds to the theoretical prediction of Eq.(||) with the 
index 9 given by Eq.(Q) and the index q determined by 
the critical value established by Fig.[| As we have seen, 
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this means q « 0.33, and, through Eq.(^), 6 w 0.4, which 
is in fact the power index of the full curve of Figjj. Note 
that the choice made to derive Eq. (^) is to some ex- 
tent arbitrary. However, the key index 9 can also be de- 
rived from 6 — 1/a, where a refers to the result of Fig|l], 
S(t) = cons t a with a ~ 2.5, coinciding with the earlier 
prediction. We also note that after the transition to the 
classical regime the decrease of ts with the increase of D 
becomes faster, as made necessary by the passage from 
the regime of Eq.(|) to that of Eq.§) . 
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FIG. 4. The saturation time ts as a function of the noise 
strength yD. The full line corresponds to the prediction of 
Eqs. (^) and ([?]) with the critical entropic index q determined 
by the numerical results illustrated in Fig.^. 

On the basis of the earlier arguments it is not a surprise 
that a condition compatible with the generalized expres- 
sion of the KS entropy is found at q ^ 1. One might won- 
der, however, why q < 1 . This can be easily explained 
by adop ting the entropic arg uments of Tsallis illustrated 
in Refs. [lq| and (l7). In fll7f the subdiffusional character 
corresponding to the entropic index q is proved by means 
of a nonlinear Fokker-Planck equation. However, there 
are good reasons to believe that this is a quite general 
property of the entropic index q 111] (see also [[P))). Thus 
we conclude that the fundamental results illustrated by 
Fig.|| and Fig.^ signal a sort of quantum-mechanically in- 
duced subdiffusion. This conclusion agrees very well with 
the well known fact that the quantum mechanical kicked 
rotor is characterized by the phenomenon of Anderson 
localization |L4]]. It is well known that the localization 
process has a rate proportional to about Ti 2 . As a conse- 
quence we expect that the entropic index q as a function 
of h tends to 1 with decreasing Ti, in agreement with 
the numerical observation that the values h = 0.10,0.05 
and 0.01 yield q = 0.28, 0.30 and 0.33, respectively. We 
thus reach the interesting conclusion that even the pro- 
cess of Anderson localization is compatible with a ther- 
modynamical treatment, provided that a proper thermo- 
dynamical indicator is adopted: This is the Tsallis en- 
tropy. Furthermore the adoption of the SE method of 
calculation makes it possible to show that if the Liou- 



ville density of Eq.(g) is coarsened by means of a weak 
stochastic force and the proper choice of the entropic in- 
dex q is made, then the time derivative of S q (t) becomes 
constant, as prescribed by the generalization of the KS 
entropy proposed by the authors of Ref . p] . 
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result:lim A:z: (o)_o = [1 + -V(l - <7)]~- Note that 

this prescription for q > 1 implies a divergence at finite 
time, and, consequently a departure of initially very close 
trajectories from one another faster that an exponential 
departure. The case q < 1 yields a power law depen- 
dence on time, which can be interpreted as slower than 
the exponential departure. It is well known that the case 
q — 1 results in ordinary Brownian diffusion. In Ref. |lij] 
it was shown that the case q > 1 leads to superdiffusion. 
Consequently, the case q < 1, with a departure of the 
trajectories slower than exponential, ought to result, in 
general, in subdiffusion, regardless whether the nonlinear 
diffusion equation of Ref. jl7j applies or not. 
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